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Abstract. We propose a model of parameter learning for signal transduction, where the objective function 
is defined by signal transmission efficiency. This is a novel approach compared to the usual technique of 
adjusting parameters only on the basis of experimental data. We apply this to learn kinetic rates as a 
form of evolutionary learning, and find parameters which satisfy the objective. These may be intersected 
with parameter sets generated from experimental time-series data to further constrain a signal transduction 
model. The resulting model is self-regulating, i.e. perturbations in protein concentrations or changes in 
extracellular signaling will automatically lead to adaptation. We systematically perturb protein concentra- 
tions and observe the response of the system. We find fits with common observations of compensatory or 
co-regulation of protein expression levels in cellular systems (e.g. PDE3, AC5). In a novel experiment, we 
alter the distribution of extracellular signaling, and observe adaptation based on optimizing signal trans- 
mission. Self-regulating systems may be predictive of unwanted drug interference effects, since they aim to 
mimic complex cellular adaptation in a unified way. 

1 A transfer function approach to signal transduction 

Signal transduction systems are usually modeled as networks of biochemical kinetic equations im- 
plemented as continuous-time dynamical models using differential equations [1J. If we regard a 
subset of species as inputs, and make sure that the system always converges to equilibrium values 
by using weakly reversible equations, we may transform these models into a set of matrices ful- 
filling the role of input-output transfer functions, i.e. a mapping from rectangular input signals 
of sufficient duration to steady-state concentrations for all target species [7]. We may then ask 
the question of optimal or efficient signal transmission for a variety of input scenarios, focusing 
on selected target species as relevant outputs or actuators of the system. An experimental anal- 
ysis has shown "fold-change" responses to changes in input [5]. A theoretical analysis has shown 
that input-output transfer functions follow the shape of hyperbolic (saturating) curves, which are 
equivalent to sigmoids for logarithmically scaled input [Tj- We hypothesize that an efficient signal 
transmission would maximize the response coefficient Rc,s (the response of species C to input S) 
defined as 
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with concentration change of target C and input S from baseline (t=0) to signal time t. (We have 
used this both with and without log, and report on logarithmically scaled input here, since the 
results seemed to be more biologically relevant.) We may also optimize for the slope s of the sigmoid 
at half-maximum concentration. This is equivalent to maximizing Rc,s, provided that the input 
signal remains entirely between the upper and lower boundaries of the sigmoid. However, shifting 
the sigmoid function to the left or to the right is even more important than slope in experimental 
settings. E.g., Fig. Q] A,B show the effect of a KO for the protein RGS in an experimental assay in 
yeast [9], and in our model. In addition to optimizing for signal transmission efficiency, we discuss 
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Figure 1: A. Experimental Response to RGS KO in yeast B. Response to RGS KO in the current model 
C. Reaction System with Selected Input(red) and Output(blue) 

optimizing for reaction time (delay to steady-state) as well. It is possible to add s as a separate 
objective. 

2 Delay vs. Efficiency Trade-off 

Fig. [Tp shows the example system, a submembrane compartment with a GPCR (G-protein coupled 
receptor) pathway from a mouse embryonic fibroblast, with ISO as input (extracellular ligand to 
beta-adrenergic receptors) and the phosphorylation of a membrane protein VASP as output. This 
model was implemented as an ODE model with 23 reactions and 10 total (initial) concentrations, 
and I/O transfer functions were derived as in [7J. The parameters were adapted to experimental 
biological data (not shown, [2]). In Fig. [2]A the original ISO/pVASP transfer function is shown, 
which we use here for further optimization. It is clear that a complex signaling system may have 
several inputs, and a large number of outputs or target proteins. It is easy to extend the present 
discussion to optimize for multiple outputs in parallel, and a system could also be optimized for 
a number of I/O functions. In that case other measures, such as mutual information, may be 
employed as objectives. The basic principle will be explained here for a single I/O pathway. The 
system was trained for a signal distribution between lOnM and 1 /iM by adjusting only the kinetic 
rate parameters. Fig. [2]A shows the I/O transfer function before and after adjustment, maximizing 
for Rc,s and/or delay. We see that optimizing for the response coefficient alone shifts the function 
to the right to better cover the input range. Optimizing for the delay alone shifts it to the left, 
speeding up signals in the lower range (which are slower). This corresponds to the theoretical results 
in [7J. Both objectives together (equally weighted) are close to the original, biologically validated 
curve. The optimization method used is a simplex algorithm [6]. These results are simple, intuitive 
and encourage continuing to explore the basic idea. 

3 Optimization and Evolutionary Learning 

In principle, we may use any parameters in a system, concentrations or kinetic rates, to maximize 
signal transmission. For the evolutionary learning which produces a cellular system, it is the kinetic 
rates which are subject to evolution of protein structure and interactions, while concentrations may 
be regulated adaptively in each cell. Nonetheless there are expected concentration ranges for healthy 
cells, specific by cell type, which are the basis for kinetic rate optimization. During evolution, new 
protein subtypes develop with a different set of interactions. This corresponds to an adjustment of 
the available set of reactions, a type of structural learning to overcome the bottlenecks that are a 
result of tightly specified molecular kinetics. 

In other words, we assume that kinetic rate learning operates on evolutionary time-scales, 
and that biologically attested signal transduction pathways contain sets of reactions and reaction 
rates which are optimal in terms of signal transmission efficiency. Our experiments show that 
kinetic rate solutions are not unique, yet there are clear tendencies on how to improve fitness of 
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Figure 2: A. Transfer Functions after Optimization of Kinetic Rates by R p vasp,isOi delay B. Relative 
Reaction Parameters during Optimization C. Binding Coefficients for High (red) and Low (blue) Fitness 

a specific I/O signal transduction pathway (Fig. [2)3). The initial activation of PKA by cAMP to 
form PKA2cAMP, the binding of ligand to receptor, the action of the GTPase RGS, the activity 
of PDE are all highly relevant, while phosphorylation of the receptor or of the cAMP synthesis 
enzyme AC are less relevant. This is evidenced by the degree of deviation of a reaction rate from 
the original as fitness improves (from left to right). (These are results from different experiments 
with different starting configurations, which were sorted according to R value. The delay values are 
not shown.) Again, this is encouraging, since it corresponds to accumulated biological observations 
on the relative importance of e.g. GTPases or PDE versus e.g. feedback phosphorylation by PKA. 

Fig. [2D shows the distribution of binding parameters according to fitness value (R, delay) of the 
signal transduction pathway. We see a distinct clustering of binding coefficients, pointing to optimal 
ranges for these parameters. It is well-known that matching models to experimental time-series 
data yields sometimes large ranges of suitable parameters [8]. 

This tool can now be used in two different ways (a) in order to optimize signal transduction using 
experimental boundaries on the rate parameters as constraints (b) in order to measure the signal 
transmission efficiency in a biologically validated system [3] and compare it to what is theoretically 
possible given the set of equations ('motifs'). 

4 Co-regulation as an Adaptive Response 

Most significantly, we now have a tool to investigate perturbations on concentrations. This is par- 
ticular interesting in the following sense: usually, when we use a computational model to investigate 
the effect of perturbation, we only study the effects. Now we can make the system itself adjust to 
the perturbation by utilizing optimization. We have created a self-regulating system, which uses 
objective functions to adapt to perturbation. This basic idea is extremely powerful, and could be 
used with different kinds of constraints on parameters, reaction times, etc. and with different I/O 
submodules of a larger system (such as cellular compartments). Here we are looking at a simple 
case (Fig. [3K): we reduce the level of a given protein (such as a kinase) to 10% of its usual con- 
centration. As a result, all concentrations in the system are allowed to adjust until optimality of 
signal transmission and delay is reinstated. This process mimics co-regulation of protein expression 
in cellular systems, which is important in disease progression and often detrimental to targeted 
interventions. The model shows strong responses to reduction of the output protein VASP, and 
the phosphatase (PP1), which prevents pVASP (phosphorylated VASP) accumulation. Lowering 
PP1 lowers most other proteins, to re-establish sensitivity, lowering VASP increases most other 
concentrations to re-establish a clear signal. PDE4B and AC6 are co-regulated in the sense that 
both are reduced together. This is interesting because it concurs with the biological facts. AC6 
synthesizes cAMP and PDE4B degrades it. I.e. the basic goal of the system is to provide for 
optimal signal transmission, here by homeostasis for the messenger molecule cAMP. Mostly it is 
only a small set of proteins that respond with altered concentrations to a single perturbation. This 
result confirms a view of signal transduction as highly modular [TJ. Concentration changes may be 
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Figure 3: A. Response to Reduction to 10% of Original Concentration B. Transfer Function Response to 
Shifts in Input Signal Range C. Concentrations Changes in Response to Extracellular Signal Shift 

caused by genetic up- or down- regulation, secretion and re-uptake, increased degradation, RNA 
interference, etc. and are therefore not easy to model from the standpoint of mechanistic biological 
modeling, which would need modules for all biologically attested processes. A unified perspective 
by a set of constraints and a set of objectives, such as has been envisaged here, may lead to better 
predictive results and may also be used as a guiding principle in constructing mechanistic models. 

5 Optimal signal transmission depends on the signal 

From the standpoint of disease modeling, an unusual protein concentration may be an adaptive 
response where a still functioning cell in a dysfunctional external signaling environment struggles 
to keep signal transmission efficient to support cellular function. In such a case, targeting this 
protein by pharmacological intervention will lead to co-regulation on other proteins. In general, 
as well as in real biological signaling systems, it is the localization of the input signaling range 
and the distribution of signals that the system transmits which are important for optimization. If 
we select the signal set that we optimize for in the right way, the input range will move towards 
the loglinear range, i.e. between the lower and upper input boundaries (Fig. [3j3) . As a result, a 
number of internal concentrations will change. Fig. [3]C shows the relative concentration changes 
that result from a shift in input range. Interestingly, the low shift requires mostly a reduction 
in RGS. The high shift (which was also less effective) affects most concentrations, notably PKA. 
The results of this process may now be compared with biological evidence, e.g. for addiction, 
where this kind of "sensitization" is well-attested, but also for cancer where intercellular signaling 
is affected. Gene expression data may then be mined not only for evidence of the mechanics of 
genetic regulatory pathways but also for evidence of shifts in extracellular signaling, which cause 
altered protein expression. 

6 Summary 

The idea to look for parametric optimization as the basis of realistic cellular properties has been 
applied with success to metabolic fluxes [5J. In that case, optimization of growth is usually regarded 
as the single objective function. Here we use another objective function, signal transmission effi- 
ciency, to study the adaptive response of a signaling system to perturbation in concentrations. This 
equals maximization of concentration change in a target species in response to input signal. Opti- 
mization of growth and optimization of signal transmission are therefore related. We also optimize 
for the speed of the transmission which counteracts signal efficiency due to the basic properties of 
kinetic equations [7J. We use optimization in a two-step process: (a) in evolutionary learning, this 
is used to find kinetic rates for estimated concentrations (b) in cellular adaptation, this is used to 
re-calibrate the model in response to perturbations in concentrations or changes in extracellular 
signaling. 
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